Phases of Holographic Superconductors in an External Magnetic Field 
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We study a 2+1 dimensional model of superconductors using a 3+1 dimensional gravitational dual 
theory of a black hole coupled to a scalar field, with negative cosmological constant. In the presence 
of finite temperature T and a background magnetic field B, we use numerical and analytic techniques 
to solve the full Maxwell-scalar equations of motion in the background geometry, finding non-trivial 
localized solutions that correspond to condensate droplets, and to vortices. The properties of these 
solutions enable us to deduce several key features of the {B, T) phase diagram. 
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I. INTRODUCTION 

For slightly over a decade, there has been a grow- 
ing branch of string theory research that extracts im- 
portant physics of strongly coupled systems by comput- 
ing with a weakly coupled gravitational "dual" system. 
These "gauge/gravity" dualities — termed thus since the 
strongly coupled system is a gauge field theory or gener- 
alization thereof — are said to be holographic in nature, 
since the dual gravitational system crucially has at least 
one extra dimension, and much of the field theory's prop- 
erties can be extracted by working on the boundary of 
the gravitating spacetime. The best understood version 
of this is the AdS/CFT correspondence [TJ [51 13] between 
superstring theory on AdSs x S'^(4+l dimensional anti- 
de-Sitter times a five-sphere) and 3+1 dimensional A/'=4 
supersymmetric SU{N) Yang-Mills theory. The theory 
can be reduced on the to a gravity theory coupled to 
a family of fields in the AdSs. The boundary of AdSs 
is a copy of 3+1 dimensional Minkowski space, and this 
is the spacetime where the dual Yang-Mills theory re- 
sides. Its 't Hooft coupling, X—ingsN—gy^N, is large 
in the regime where the dual gravity model is reliable: the 
large N limit and weak string coupling gg. The asymp- 
totic values of the AdSs fields supply information about 
operators in the Yang-Mills theory. Ref.jlj reviews of 
much of this technology, with bibliography. This type of 
duality has supplied a wealth of information about vari- 
ous strongly coupled systems that can be supplied with 
a gravitational (or fully stringy) dual along these lines, 
ft has several potentially important applications, rang- 
ing from models of the behaviour of the strongly coupled 
dynamics of quark-gluon plasmas in nuclear physics, to 
models relevant to condensed matter physics, all of con- 
siderable experimental interest. (See refs. [Sj [6l [7] for re- 
views and bibliography.) 

A holographic model of some of the key phenomcno- 
logical attributes of superconductivity in 2+1 dimensions 
was proposed in ref.[8j. The dual is a simple model of 



gravity in four dimensions (AdS4) coupled to a C/(l) 
gauge field and a minimally coupled charged complex 
scalar ^ with potential V{\^\) = ~2\^\'^/L'^, where the 
cosmological constant defines the scale L via A = 
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where k\ — SttGn is the gravitational coupling and our 
signature is ( — h ++). We will use coordinates (t, z, r, 0) 
for much of our discussion, with t time, (r, 0) forming a 
plane, and z a "radial" coordinate for our asymptotically 
AdS4 spacetimes such that z = is the boundary at in- 
finity. Note that the mass of the scalar to|, = — 
is above the Brcitenlohner-Freedman stability bound [!Jj 
m|p = — 9/4L^ for scalars in AdS4. A black hole (which 
is planar, its horizon is an (r, cj)) plane at some finite 
z — Zh), with Hawking temperature T and mass per unit 
horizon area e = M /V , corresponds to the dual 2+1 
dimensional system at temperature T and with energy 
density e. The asymptotic value of ^' on the boundary 
sets the vacuum expectation value (vev) of a charged op- 
erator O in the system, which is the order parameter. 
For T > Tc, the scalar (and hence (O)) is zero, and for 
T < Tc it is non-zero. In the gravitational theory, the 
high T phase is simply the charged black hole (Reissner- 
Nordstrom, or AdS-RN) with vanishing. Notice that 
the mass of the scalar is set not just by T^d^l/I) but by 
the black hole's gauge field A = Atdt. {At does not 
give an electric field in the dual theory on (r, 0) , but de- 
fines instead a U{1) charge density [TP], p. See below.) 
T also depends on p. In fact to|, decreases with T until 
at Tc it goes below r7i|p, becoming tachyonic. The the- 
ory seeks a new solution, in which the black hole is no 
longer AdS-RN, but one that has a non-trivial profile 
for (i.e., it has "scalar hair" — for a discussion of vio- 
lations of non-hair theorems in this context, see ref. jlli). 
The U{1) is broken by (O). These solutions can be found 
by solving equations of motion in certain limits (explored 
below) , and the transport properties of the low tempera- 
ture phase were examined in refs.[8| [T2]. using linear re- 
sponse theory, with the result that the DC conductivity 
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diverges in a way consistent with expectations that the 
phase is superconducting. (Strictly speaking, the C/(l) 
that was broken by (O) is global on the boundary, but it 
can be gauged in a number of ways without affecting the 
conclusions. See e.g. ref.[12j.) 

It is clearly of interest to study this system further, 
since it could well open the door to a whole new phe- 
nomenology of superconductivity in a wide range of phys- 
ical systems. We report here on our study of the sys- 
tem in an external magnetic field, continuing the work 
we began in ref. [13]. The magnetic field B (which fills 
the two dimensions of the superconducting theory) , also 
contributes to m|,, via its square, but contributes with 
opposite sign to the electric contribution of the back- 
ground. It therefore lowers the temperature at which 
to|, falls below m^p, triggering the phase transition. On 
these grounds alone one then expects a critical line in the 
{B, T) plane of a form like that of a full curve from fig. |3] 
but it is important to determine exactly what physics lies 
on either side of the line. Generically, for non-zero B, it 
is inconsistent to have non-trivial spatially independent 
solutions on the boundary, and we study two classes of 
localized solutions in some detail. The first is a "droplet" 
solution, the prototype of which was found in our earlier 
work as a strip in 2D (straightforwardly generalized 
to circular symmetry in ref . fl2' ) , and the second is a vor- 
tex solution, with integer winding number ^ G -E, which 
is entirely new. We obtain these as full solutions of the 
Maxwell-scalar sector in a limit, and their properties al- 
low us to determine key features of the {B,T) phase di- 
agram (correcting statements made in ref. |12j'). 



choose to study the case of i=l, for brevity. ^ sets a 
chemical potential, and p sets a charge density for the 
U{1) symmetry that we are studying here. To find a 
full solution of the equations, we solve the equations of 
motion numerically, using a shooting method starting at 
the horizon (2:=!). There we set p{l) to a constant, 
and for regularity of the gauge field, we have At{l)=0. 
To look for solutions we tune dzAt{l) such that at the 
other boundary, p(0) has the required Neumann condi- 
tion p2=0. From that we read off pi=y/2K4{Oi) /{La), 
and p. Since the only scale in the theory besides T is set 
by the charge density p, the value of Tc is given in terms 
of p. This can be determined by noting when a non-zero 
(Oi) develops, and there we find Tc—0.226ay/p. Fig. 1 
shows the result for the vev as a function of temperature, 
showing the low and high temperature phases separated 
by a second order phase transition at Tc below which the 
U{1) is spontaneously broken. 

Tc 
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FIG. 1: The scalar p ~ (d) vs. T. Tc = 0.226q^. 



II. SPATIALLY INDEPENDENT SOLUTIONS 

This section reviews the spatially independent conden- 
sate solution of ref. [8J that corresponds to the supercon- 
ducting (symmetry breaking) phase below some T^. It 
is found in a certain decoupling limit. Under a redefi- 
nition An/g, ^'/g, the Maxwell-scalar part 
of the action M gets a prefactor of l/g^, and the g dis- 
appears from the A'i' coupling. This means that in the 
limit 5 — > 00, this sector decouples from gravity. We 
can therefore take as background black hole the (planar) 
Schwarzschild solution: 

ds' = ^ {~f{z)de + dr^ + r'dd^') + -^^dz' , 

(2) 

with f{z)=l — z^ . It has an horizon at z—1, temper- 
ature T=3a/47r, and mass density £—L^a^/K\. Let us 
write our complex scalar as ^—peyi\i{i9)/^/2L and write 
At=aAt. The equations of motion allow an ansatz: 
6'=const., p—p{z), At—At{z), and A^=0. Near the 
boundary z—0 we have, for constants pi, p2, p, p'. 

p^piz + p2Z^, At^p~pz. (3) 

Pi (i=l, 2) sets the vev of a A=i operator Oi [Tl|. Only 
one of these vevs can be non-zero at a time, and we will 



III. SPATIALLY DEPENDENT SOLUTIONS 

We can seek non-trivial solutions with 6=(^ + ^0, and 

p = p{f, z),At^ At{f, z),A^^ A^{r, z) , (4) 

where r—ar and ((., 1^) are constants, with ^ integer. Reg- 
ularity of the equations of motion require that near f=0, 
we must have p^f^ . This feeds into the behaviour of all 
the other fields near f=0, and determining this near the 
horizon is important in order to seed the numerical search 
in a manner that will carefully find solutions. We report 
on this in a longer publication[15 . The solutions that 
we find using our analysis break into two broad classes. 
There are those for which p — )■ for f ^ 00 which are 
called droplets. The prototype localized solution of this 
type was found in our earlier workjl3j. and also stud- 
ied in ref. [12], but in a "probe" limit of small fields and 
finite g. Here we have the full solutions in the large g 
decoupling limit. Those solutions for which p — > const, 
for f ^ 00 we call vortex solutions (^ G 2+). They be- 
have exactly as expected of vortices in this context. Note 
that ^ defines a non-trivial topological winding number. 
A^ becomes constant at infinity so gauge symmetry can- 
not be used to unwind 0. Gauge symmetry is unbroken 
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at infinity for the droplets, so ^ is not a winding number 
for them. 



A. The Droplet 

A sample droplet solution is presented in fig. [2] (for p 
and Bz). The current density and magnetic field are 
read off via the z ^ asymptotic — > + J^(r)z, 
with Bz=dfa^/r. The constant 7 appears in the small r 



and has as solution: i?(f)=i?o(l) exp (— 7f^/4) . We use 

this as the seed to solve the coupled equations for At and 
R for the full z dependence in this limit, seeking the ap- 
propriate boundary conditions at z=0. This extracts, as 
in section |TT] the T of the solution for our magnetic field. 
This procedure sets the critical magnetic field at which 
the droplet solutions first form for a given T/Tc- Be- 
low that B the droplets simply vanish since in this probe 
limit they are already at zero size. Our result is the top 
curve in fig. [3] Large fields will begin to back-react on 
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FIG. 2: A droplet solution for 7 = 0.5 at T/Tc ~ 0.84. 

expansion of the function R defined hy p — zr^R, where 
^ ^ i?o (1 — 7r^/4 + • • • ) at the horizon. We see that 
the magnetic field fills the plane, asymptoting to a con- 
stant value at r —^ 00, and approaching a non-zero value 
in the core. We studied several solutions for a range of 
T and B = B^, noticing that for high temperatures B 
is decreased somewhat in the core as compared to the 
asymptotic value, while for low temperatures B is en- 
hanced there. 

It is important to determine exactly where in the 
{B,T) plane these solutions can appear. We find our 
full solutions at a variety of values of B, but none below 
certain values of B, which is suggestive. To determine if 
there is a minimum value of B to form droplets, we can 
use the probe limit, where p is small. Here we can take 
A,f, = 7f^/2 (which will determine the magnetic field in 
this limit) and have At — At{z). The equation of motion 
for R on the horizon reduces to: 



diR + -dfR - -P^j^R + jR 
r 4 



= 



(5) 
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FIG. 3: Scaled {B, T) diagram connecting to decoupling limit. 

the geometry, and our analysis breaks down, as can be 
seen in our curve for low T/Tc- To proceed further, we 
can carry out the probe computation again, but not in 
the g ^ 00 limit. Since the Maxwell sector is not decou- 
pled from the geometry for finite g we use as background 
a dyonic black hole. This is the analysis of our earlier 
work [13]. The metric is of the form given in eqn. Q, 
but with a Maxwell field F = 2ha'^rdr t\d(j) + 2qadz A dt 
and f{z) = (1 - z) (z^ + ^ + 1 - (/i2 + ^3)^ Nq^^ y 

and p are determined by the background (we do not dis- 
play them here). We note that the equations of mo- 
tion have a separable solution p ~ zZ{z)R{r), and R{f) 
satisfies an equation of the same form as ([5|, but with 
7 2gh, giving the Gaussian profile. We can then 
solve for Z(z) numerically using the same shooting tech- 
niques as before. This gives a complete {B, T) curve for a 
given g. In fact, these curves connect to the g ^ 00 probe 
computation of above, as can be established by careful 
comparison of the different definitions of the tempera- 
tures and charge densities in the two limits. We carry 
this out in our longer paper [TS]. and simply show here 
how the two limits connect by rescaling our curves (with 
the appropriate factor of g) and superposing. See fig. [3| 



B. The Vortex 

We display the ^ = 1 vortex solution in fig. |4] (for p 
and Bz). As noted before, there is one unit of winding, 
and the scalar field runs to a constant at infinity. So does 
the charge density, and so we are able to read off the value 
of the temperature for these solutions. The magnetic field 
drops to zero at large f, and A^ becomes constant there. 
That constant is ^ in general, as is consistent with the 
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total magnetic flux being quantized to 27r^. The current 
density Jcj,{f) (not displayed for brevity) is zero asymp- 
totically and peaks in a ring around the core, supporting 
the magnetic field, as expected for a vortex. As for the 



P(?) 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



0.30t"\ 
0.25 \ 
0.20 \ 
0.15--N^\ 
0.10 \.» 
0.05 



T 

— = 0.92 
— Tc 

T 

— = 0.77 
.... J 



6 8 10 12 14 
(a) Scalar Field 



V 



T 






= 0.92 


T, 




T 






= 0.77 


t;- 





2 4 6 8 10 12 14 
(b)Magnetic Field 

FIG. 4: The vortex solution for ^ = 1. 



droplets, the question arises as to where in the {B,T) 
plane vortices appear. As mentioned, we determine T/Tc 
from the large f values of p and p\z=o — V^K4{Oi) /{La), 
while the magnetic field is localized in the core. For an 
applied B at some T/Tc, we expect that the system cre- 
ates vortices to trap the flux in their cores, forming a 
lattice of them that increases in density with B. Vortices 
presumably repel each other, and so such a lattice will 
cost energy. Therefore at some critical Bc{T) the system 
will seek a lower energy phase, possibly returning to the 
normal phase. We do not know if this critical line co- 
incides with the one in fig. [3] where the droplets begin, 
although it is the simplest possibility. 



IV. CONCLUSION 

We have studied the prototype holographic supercon- 
ductor of ref.[H] in the presence of magnetic field, con- 
tinuing our earlier work|13j. examining two important 
classes of solution, the droplet and the vortex. The 
latter is completely new, and this is the first time the 
droplet has been fully constructed. Our analysis in var- 
ious connected limits shows where these solutions can 
exist in the {B,T) plane. There is a critical line below 
which droplets are not found, while vortices can be found 
there. Our interpretation is that this region is the super- 
conducting phase, and that for non-zero B, the vortices 
develop, trapping the magnetic flux into fllaments, as is 
familiar in type II superconductors. Above the critical 
line, the system leaves the superconducting phase, and 
either forms droplets of condensate or simply reverts to 
the normal phase (dual to a dyonic black hole with zero 
scalar everywhere, which may well yield lower action than 
the droplets if we had back-reacting solutions to work 
with). We have performed a stability analysis of these 
solutions [T5|. Both classes are stable for the examples 
and modes that we study (this is expected for the vor- 
tices given their conserved winding). Finally, note that 
we disagree with the suggestion made in ref. [T^]. The 
authors find the critical line, but state (similarly to our 
ref.[T3]) that the droplets exist below the line, and are 
superconducting. As they did not have the full droplet 
solutions, nor the vortex solutions, their analyses are not 
sufficient to make these determinations. Our work here 
gives a much stronger picture of the phase diagram. 
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